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Preface 


The  purpose  of  this  thesis  is  to  generate  critical 
values  for  the  modified  Anderson-Darling  and  Cramer-von  Mises 
statistics.  These  critical  values  are  used  for  testing 
whether  a  set  of  observations  follows  a  Weibull  distribution 
when  the  scale  and  location  parameters  are  unspecified  and 
are  estimated  from  the  sample.  An  extensive  power  study  is 
made  to  compare  the  power  of  the  Anderson-Darling,  Cramer- 
von  Mises,  Kolmogorov-Smirnov,  and  the  Chi-Square  goodness- 
of-fit  tests. 
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Abstract 

The  Anderson-Darling  and  Cramer-von  Mises  critical 
values  are  generated  for  the  three-parameter  Weibull  distri¬ 
bution.  The  critical  values  are  used  for  testing  whether  a 
set  of  observations  follows  a  Weibull  distribution  when  the 
scale  and  location  parameters  are  unspecified  and  are  esti¬ 
mated  from  the  sample.  A  Monte  Carlo  simulation,  with  5000 
repetitions,  is  used  to  generate  critical  values  for  sample 
sizes  5(5)30  and  Weibull  shape  parameters  .5(. 5)4.0. 

A  Monte  Carlo  power  investigation  of  the  Anderson- 
Darling  and  Cramer-von  Mises  tests  is  made  using  5,  15,  and 
25  observations  from  ten  alternate  distributions.  The  power 
of  the  two  tests  are  compared  to  the  Kolmogorov-Smirnov  and 
the  Chi-Square  tests.  The  power  of  all  the  tests  are  low 
with  a  sample  size  of  five.  When  the  hypothesized  distribu¬ 
tion  is  the  Weibull  with  shape  equal  1.0,  the  power  of  the 
tests  in  decreasing  order  are:  Cramer-von  Mises,  Anderson- 
Darling,  Kolmogorov-Smirnov,  and  Chi-Square.  When  the 
hypothesized  distribution  is  the  Weibull  with  shape  equal 
3.5,  the  power  of  the  tests  are  the  Anderson-Darling,  followed 
by  the  Cramer-von  Mises,  Kolmogorov-Smirnov,  and  the  Chi-Square 
test. 

A  relationship  between  the  Anderson-Darling  and 
Cramer-von  Mises  critical  values  with  the  Weibull  shape 
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A  MODIFIED  CRAMER-VON  MISES  AND 
ANDERSON -DARLING  TEST  FOR  THE 
WEIBULL  DISTRIBUTION  WITH 
UNKNOWN  LOCATION  AND 
SCALE  PARAMETERS 

I .  Introduction 

Reliability  theory  and  life  testing  of  equipment  is 
of  great  importance  to  the  U.S.  Air  Force.  The  mean  time  to 
failure  and  the  failure  rate  ot  equipment  and  systems  are  an 
important  input  into  the  decision  making  process. 

Through  experimentation  and  observation,  time  to 
failure  data  can  be  collected.  The  data  can  be  compared  to 
a  theoretical  probability  distribution.  The  test  to  deter¬ 
mine  if  the  hypothesized  distribution  "fits"  the  data  in  the 
sample  is  known  as  the  "goodness-of-fit  test."  There  are 
several  classical  goodness-of-fit  tests,  such  as:  the  Chi- 

Square  test  (x  ) ,  the  Kolmogorov-Smirnov  test  (K-S) ,  the 

2 

Anderson-Darling  test  (A  ) ,  and  the  Cramer-von  Mises  test 
2 

(W  ) .  If  such  tests  show  a  good  "fit,"  the  hypothesized  dis¬ 
tribution  can  then  be  used  in  simulation  modeling  techniques 
to  predict  the  failure  rates  of  Air  Force  systems  or  system 
components . 


1 


Background 


The  general  goodness -of-fit  tests  are  valid  to  test 
whether  a  set  of  observations  come  from  a  completely  speci¬ 
fied  distribution.  The  tests  have  been  modified  when  the 
parameters  of  the  distribution  are  unknown  and  are  estimated 
from  the  data  sample.  The  K-S  test  has  been  modified  by 
H.W.  Lilliefors  so  the  test  can  be  used  with  the  normal  (16) 
and  the  exponential  (17)  distributions  when  the  parameters 
are  unknown  and  are  estimated  from  the  data.  R.  Cortes 
further  extended  the  K-S  test  so  it  can  be  used  with  the 
Gamma  and  Weibull  distributions  when  the  scale  and  location 

parameters  are  unknown  (4).  In  1969,  Green  and  Hegezy  (10) 

2  2 

used  the  K-S,  W  and  the  A  tests  to  generate  rejection  tables 
for  the  goodness-of-fit  for  the  uniform,  normal,  Laplace, 
exponential,  and  the  Cauchy  distribution  with  unknown  para¬ 
meters.  Mann,  Scheuer,  and  Fertig  developed  in  1973  a  new 
goodness-of-fit  test  for  the  two-parameter  Weibull  distribu¬ 
tion  with  unknown  parameters.  They  called  the  new  statistic 

the  L  and  S  statistic  (19).  In  1979,  Littell,  McClave,  and 

2  2 

Of fen  (18)  used  the  K-S,  W  ,  and  the  A  tests  to  generate  the 
rejection  regions  for  the  two-parameter  Weibull  distribution 
when  the  shape  and  scale  parameters  are  unspecified. 

Empirical  Distribution 
Function 

The  empirical  distribution  function  (EDF)  statistics 
are  a  class  of  statistics  based  on  a  comparison  between  the 
cumulative  distribution  function  F(x)  and  the  empirical 


2 


distribution  function  Sn(x) .  The  EDF  of  a  random  sample  of 
size  n  is  the  proportion  of  the  sample  values  which  do  not 
exceed  the  number  x: 

Sn(x)  *  number  of  sample  values  1  x  ^ 

Thus,  Sn(x)  is  a  step  function  with  1/n  jumps  at  each  order 
statistic  of  the  sample  (9:73).  Letting  x^ ,x(2) * ’ * *x(n) 
be  order  statistics,  Sn(x)  is  defined  by  Eq  (2). 


Sn(x) 


'0  ,  x  <  x(1) 

S  i/n  *  x(i)  <  *  <  x(i+1) 
-1  ’  x  >  x(n) 


1*1,2 


(n-1) 


(2) 


Since  Sn(x)  gives  the  proportion  of  a  random  sample  less  than 
x,  it  is  reasonable  to  expect  it  to  give  a  good  estimate  of 
the  hypothesized  cumulative  distribution  function  F(x) .  In 
fact,  goodness-of- fit  tests  based  on  EDF  statistics  measure 
the  discrepancy  between  Sn(x)  and  F(x). 

The  EDF  statistics  used  in  this  thesis  are  the 
Anderson-Darling  statistic  and  the  Cramer- von  Mises  statistic . 
They  are  compared  with  perhaps  the  best  known  EDF  statistic, 
the  Kolmogorov- Smirnov  goodness-of-fit  test.  The  K-S  statis¬ 
tic  measures  the  maximum  vertical  discrepancy  between  the 
cumulative  distribution  function  and  the  empirical  distribu¬ 
tion  function  (23:2). 

Cramer-von  Mises  Goodness- of -Fit  Test  Statistic.  The 
Cramer-von  Mises  statistic  is  based  on  the  squared  integral  of 
the  difference  between  the  cumulative  distribution  tested  and 


3 


the  EDF  (23:2)  : 


W*  =  £  [Sn(x)  -  F(x) ] 2  ( x) dx  (3) 


where  ijj(x)  gives  a  weighting  to  the  squared  difference.  The 

Cramer- von  Mises  statistic  is  figured  with  if»(x)  *  1.  The 

2 

computational  formula  to  calculate W  is  given  by  Eq  (4). 

Letting  x(i)  »x( 2)  * '  *  *x(n)  t*ie  n  order  statistics 
and  letting  IK  *  F(x^) ,  the  cumulative  distribution  function 
then, 


W 


2 


n 


T2n  +  if1[Ui 


2i-l,2 
2n  1 


(4) 


2 

If  W  is  too  large,  the  hypothesized  distribution  is  rejected 
(1:765) . 

Anderson-Darling  Goodness-of-Fit  Test  Statistic.  The 

Anderson-Darling  goodness -of -fit  statistic  is  also  based  on 

the  squared  integral  of  the  difference  between  F(x)  and  Sn(x) . 
7 

The  A  statistic  is  given  by  Eq  (3]  with 


<Kx)  *  (CF(x))(l-F(x))]'1  (5) 

This  weight  function  has  the  effect  of  giving  greater  import¬ 
ance  to  observations  in  the  tails  (23:2).  The  computational 

2 

formula  to  calculate  A  is  given  by  the  next  equation  (16) : 

A2  -  *n  -  i{Ji(2i-l)[ln(Ui)+ln(l-Un+1..)]}  (6) 

2 

If  A  is  too  large,  then  the  hypothesis  is  rejected  (23:4). 


4 


Unknown  Parameters.  Because  of  the  probability  inte¬ 


gral  transformation,  the  values  of  a  completely  specified 
cumulative  distribution  are  ordered  values  from  a  uniform 
distribution  over  the  interval  from  zero  to  one.  EDF  goodness 
of- fit  tests  test  the  hypothesis  that  a  sample  has  been  draivn 
from  a  fully  specified  continuous  cumulative  distribution. 
Thus,  EDF  statistics  are  a  function  of  ordered  uniform  random 
variables  (23:4).  In  general,  if  the  cumulative  distribution 
function  is  not  completely  specified,  the  distribution  of  the 
EDF  statistics  will  depend  on  the  sample  size  n  and  the  values 
of  the  unknown  parameters  (23:4).  However,  David  and  Johnson 
showed,  in  1948,  that  the  distribution  of  any  EDF  statistic 
is  simplified  when  the  unknown  parameters  are  the  location 
and  scale.  They  showed  the  distribution  of  any  test  statis¬ 
tic  based  on  the  cumulative  distribution  function  using  invari 
ant  estimators  for  location  and  scale  will  depend  on  the 
distribution  tested,  but  not  on  the  specific  values  of  the 

unspecified  scale  and  location  parameters  (5:182).  Thus,  the 

2  2 

distribution  of  the  A  ,  W  ,  and  K-S  statistic  for  the  Weibull 
distribution  with  unknown  scale  and  location  parameters  will 
depend  on  the  sample  size  n,  but  will  be  independent  of  the 
scale  and  location  parameters. 

Problem 

It  is  well  known  that  the  EDF  goodness-of- fit  tests 
are,  in  general,  more  powerful  tests  of  the  null  hypothesis 
than  is  the  Chi-Square  test  (24:730).  The  most  powerful 
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goodness -of- fit  tests  among  the  Anderson-Darling,  Cramer- von 
Mises,  and  the  modified  Kolmogorov-Smirnov  tests  for  the 
three-parameter  Weibull  distribution  when  the  scale  and  loca¬ 
tion  parameters  are  not  specified  is  not  known. 

Objectives 

This  thesis  has  the  following  objectives: 

1.  To  generate  the  Anderson-Darling  rejection  tables 
for  the  three-parameter  Weibull  distribution  when  the  scale 
and  location  parameters  are  not  specified. 

2.  To  generate  the  Cramer- von  Mises  rejection  tables 
for  the  three-parameter  Weibull  distribution  when  the  scale 
and  location  parameters  are  not  specified. 

3.  To  conduct  a  power  comparison  between  the  K-S, 

2  2  2 

A  ,  W  ,  and  the  x  goodness-of- f it  tests  for  the  three- 
parameter  Weibull  distribution  when  the  scale  and  location 
parameters  are  not  specified. 

4.  To  investigate  a  relationship  between  the  shape 
parameter  of  the  Weibull  distribution  and  the  Anderson-Darling 
and  Cramer- von  Mises  critical  values. 
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II.  The  Weibull  Distribution 


History  and  Application 

The  Weibull  probability  density  function  (pdf)  was 
developed  in  1939  by  a  Swedish  scientist,  Waloddi  Weibull. 

In  his  study,  he  examined  the  distribution  of  the  phenomenon 
of  rupture  in  solids.  In  a  paper  published  in  1951,  Weibull 
demonstrated  the  application  of  the  distribution  in  the  inves¬ 
tigation  of  the  yield  strength  and  fatigue  of  steel,  the  size 
distribution  of  fly  ash,  and  the  fiber  strength  of  cotton 
(26:293).  Although  the  Weibull  distribution  was  first  applied 
to  the  fatigue  of  materials,  the  distribution  is  very  flexible 
and  has  proved  very  useful  in  other  fields  of  study  as  well. 
Peto  and  Lee  used  the  Weibull  distribution  in  their  experi¬ 
ments  with  continuous  carcinogenesis  experiments  with  labora¬ 
tory  animals  (21:457),  Today,  the  Weibull  distribution  is 
best  known  for  its  application  in  the  field  of  reliability. 

The  Weibull  distribution  is  related  to  the  exponen¬ 
tial  distribution.  However,  the  Weibull  distribution  with 
the  three  parameters- -shape ,  location,  and  scale- -can  model 
a  variety  of  hazard  situations.  Where  the  exponential  has  a 
constant  failure  rate,  1/9,  the  Weibull  distribution  can  deal 
with  an  increasing,  constant,  or  decreasing  failure  rate. 

Since  many  systems  often  experience  a  burn-in  period  (decreas¬ 
ing  failure  rate)  followed  by  a  useful  life  (constant  failure 
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rate)  followed  by  a  wearout  period  (increasing  failure  rate), 
the  Weibull  distribution  can  be  very  useful  in  describing 
these  phenomena. 

The  Three-Parameter  Weibull 

The  Weibull  pdf  represents  the  distribution  of  x, 
such  as  the  time  of  failure  of  a  piece  of  equipment.  Let  K 
denote  the  shape  parameter,  C  the  location  parameter,  and  0 
the  scale  parameter.  The  Weibull  pdf  is: 

f  K(x-~-CIK  1  exp{-(^-)K}  ,  K,6>0  ,C<x 

eK  0 

F(x :  K,C ,  9)  =<  (7) 

0  ,  otherwise 

The  Weibull  cumulative  distribution  function  (CDF) ,  F(x)  is 
given  by: 


F(x)  -  /* 

f (x: K,C ,0 ) dx 

(8) 

»  1  ■ 

-  exp{-(^)K}  K ,  9>Q  ,  C<x 

(9) 

The  mean  and  variance  of  the  Weibull  distribution  are: 

Mean  -  C  +  0r(^i)  (10) 

Variance  *  C2[r(^)-  r2(£j^)]  (11) 

If  there  are  n  identical  equipment  components  with 
the  same  failure  distribution  F(x) ,  then  the  hazard  function 
is  the  percentage  of  the  n  components  that  will  fail  in  the 
next  time  interval  x.  The  hazard  function  is  defined  as: 
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The  hazard  function  for  the  three-parameter  Weibull  distribu¬ 
tion  is: 


h(x) 


(13) 


The  shape  parameter  K  determines  whether  the  hazard  function 
is  decreasing,  constant,  or  increasing.  A  shape  parameter  K 
less  than  unity  has  a  decreasing  hazard  rate  and  describes 
the  burn-in  period.  A  shape  parameter  K  equal  to  unity  has 
a  constant  failure  rate,  1/6,  and  describes  the  useful  life 
period.  A  shape  parameter  K  greater  than  unity  has  an  in¬ 
creasing  hazard  rate  and  describes  the  wearout  period.  This 
description  is  commonly  known  as  the  bath  tub  curve  as  shown 
in  Fig.  1. 

The  Weibull  distribution  reduces  to  the  exponential 
distribution  when  K  is  set  to  one  and  C  is  set  to  zero.  The 
Weibull  distribution  approximates  the  normal  distribution 
with  K  =  3.4  (3:45).  Figure  2  depicts  the  Weibull  with  dif¬ 
ferent  values  of  K  while  C  and  6  are  held  constant.  This 
flexibility  in  describing  different  failure  rates  makes  the 
Weibull  distribution  most  appealing. 

The  location  parameter  C  indicates  the  value  x  at 
which  failures  begin  occurring.  With  C  equal  zero,  failures 
can  occur  immediately  after  the  component  or  device  is  put 
into  operation.  A  location  parameter  greater  than  zero  indi¬ 
cates  that  there  is  a  period  of  time  which  is  failure  free. 


Deceasing  Constant  hazard  rat*  Increasing  hazard  rat* 

hazard  rata 

Operating  time  (component  age  or  life) 

Fig.  1.  Bath  Tub  Curve  (7:28) 

A  location  parameter  less  than  zero  indicates  that  the  com¬ 
ponents  are  defective  at  the  time  of  initial  operation. 
Failure  during  storage  of  the  component  would  be  modeled 
using,  a  negative  location  parameter.  In  this  thesis,  it 
is  assumed  C  is  non-negative.  In  Fig.  3,  C  is  varied  while 
the  shape  parameter  and  scale  parameters  are  held  constant. 

The  scale  parameter  0  affects  the  dispersion  of 
the  random  variable,  x,  about  its  mean  and  is  sometimes 
called  the  characteristic  life  (3:48).  Figure  4  depicts 
various  values  of  0  while  the  shape  and  location  parameters 
are  held  constant. 
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2.' 


Fig.  2.  Weibull  Distribution  with 
Different  Values  of  K;  6=1;  C=0  (2:47) 


i.oj- 


Fig.  3.  Weibull  Distribution  with 
Different  Values  of  C*.  6  =  1;  K=2  (2:48) 


o 


Fig.  4.  Weibull  Distribution  with 
Different  Values  of  8;  K*2;  C*0  (2:46) 
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Ill .  Methodology 

This  thesis  presents  a  Monte  Carlo  method  for  obtain¬ 
ing  the  critical  values  of  the  modified  Anderson-Darling  and 
the  modified  Cramer-von  Mises  goodness -of- fit  for  the  three- 
parameter  Weibull  distribution  when  the  scale  and  location 
parameters  are  not  specified.  This  chapter  discusses  the 
procedures  used  in  this  thesis.  First,  an  outline  and  flow 
chart  of  the  Monte  Carlo  method  is  presented.  This  is  followed 
by  a  detailed  description  of  the  steps  taken  in  the  Monte  Carlo 
procedure.  The  chapter  concludes  with  a  description  of  the 
power  study  and  of  the  analysis  of  the  critical  values  for 
each  goodness-of- fit  test  versus  the  Weibull  shape  parameters. 

Steps  in  the  Monte  Carlo 
Method' 

The  following  nine  steps  outline  the  Monte  Carlo 

2 

method  used  to  calculate  the  critical  values  for  the  A  and 
2 

W  goodness-of-fit  tests.  The  flow  chart  in  Fig.  5  illustrates 
the  method. 

1.  For  a  fixed  sample  size  n  and  fixed  shape  para¬ 
meter  K,  n  random  Weibull  deviates  are  generated  using  a 
computer  subroutine.  All  Weibull  deviates  are  generated  with 
location  parameter  C=2  and  scale  parameter  9=1. 

2.  The  n  random  deviates  are  ordered,  x(i)>x(2)» 

,X(n) * 
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Fig  5.  Flow  Chart 


Fig  5,  Flow  Chart  continued 


© 


Step  6 


Step  7 


Step  8 


Step  9 


3.  The  ordered  random  Weibull  deviates  are  used  to 
estimate  the  scale  and  location  parameters  by  the  method  of 
maximum  likelihood. 

4.  The  estimated  scale  and  location  parameter  and 
the  fixed  shape  parameter  are  used  to  determine  the  hypothe¬ 
sized  distribution  function  F(x). 

5.  The  Anderson-Darling  statistic  is  calculated 
using  Eq  (6) . 

6.  Steps  1  to  5  are  repeated  5000  times,  thus 
generating  5000  independent  A2  statistics. 

7.  The  5000  statistics  are  ordered.  Using  a  plot¬ 
ting  position,  described  later,  the  80th,  85th,  90th,  95th, 
and  the  99th  percentiles  are  calculated  by  linear  interpola¬ 
tion. 

8.  Repeat  steps  1  to  7  for  sample  size  n  equal  to 
5,  10,  15,  20,  25,  and  30. 

9.  Repeat  steps  1  to  8  for  shape  parameter  K  equal 
to  .5,  1.0,  1.5,  2.0,  2.5,  3.0,  3.5,  and  4.0. 

The  critical  values  for  the  Cramer-von  Mises  test  is 

2 

calculated  using  the  same  procedure  by  calculating  the  W 
statistic  in  step  5  using  Eq  (4) . 

Generation  of  the  Random  Weibull  Deviates.  The 
random  Weibull  deviates  are  obtained  on  the  Control  Data 
System  (CDC)  6600  computer  using  the  International  Mathemati¬ 
cal  and  Statistics  Library  (IMSL)  subroutine  GGWIB  (12:  Ch.  G) . 

This  subroutine  generates  pseudo- random  deviates 
by 
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X  =  [-ln(u) ] 


(14) 


1/K 

where  u  is  a  pseudo- random  deviate  from  a  uniform  (0,1)  dis¬ 
tribution.  Thus,  the  deviates  obtained  by  this  subroutine 
are  from  the  Weibull  distribution  with  location  parameter 
C  *  0,  scale  parameter  9=1,  and  the  shape  K  set  as  an  input 
parameter.  Deviates  from  the  more  general  Weibull  pdf  are 
obtained  using  the  transformation 

x*  =  0  •  x  +  C  (15) 

In  this  thesis  0=1  and  C  =  2.  The  reason  for  this  step 
will  be  explained  after  the  parameter  estimation  routine  is 
discussed. 

Ordering  of  the  Deviates.  The  random  Weibull  devi¬ 
ates  are  ordered,  smallest  to  largest,  using  the  CDC  6600 
IMSL  subroutine  VSRTA  (12:Ch.  V). 

Method  of  Maximum  Likelihood  Estimation  of  the 
Weibull  Parameters.  The  procedure  to  derive  the  maximum 

AAA 

likelihood  estimators  K,  C,  0  of  the  Weibull  parameters  K, 

C,  0  was  developed  by  Harter  and  Moore  (11:639).  This  is  a 
general  purpose  iterative  method  for  censored  or  uncensored 
samples.  The  Harter  and  Moore  parameter  estimation  routine 
maps  negative  estimates  onto  zero.  If  Weibull  deviates  are 
generated  with  location  C  »  0  and  scale  9=1,  the  estimators 

A  /V 

C  and  9  will  not  retain  the  invariant  property  because  within 
the  iterative  estimation  procedure  negative  location  estima¬ 
tors  are  mapped  onto  zero.  Thus,  Weibull  deviates  are 
generated  with  C  =  2  and  9=1,  preventing  the  generation  of 
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negative  location  estimators.  Because  of  the  invariant  pro¬ 
perty  of  the  EDF  statistics,  this  transformation  will  not 
adversely  affect  the  results. 

The  method  of  maximum  likelihood  selects  as  estimates 
those  values  of  the  parameters  that  maximize  the  probability 
or  joint  density  of  the  observed  sample  (20:302).  The  likeli¬ 
hood  function  is  defined  as  follows.  Let  x^,x2,...,xn  be 
sample  observations  taken  on  corresponding  random  variables: 

Xj ,X2 , . . . ,Xn.  Then  if  X1>X2,...Xn  are  continuous  random 
variables,  the  likelihood,  L,  is  defined  to  be  the  joint 
density  evaluated  at  »x2 , . . . ,xR  (20:303).  For  the  Weibull 
pdf,  the  likelihood  function  can  be  represented  by: 

L  ■  f(Xj,x2>..., xn :K,8 ,C)  (16) 

Since  the  deviates,  x^,  are  derived  randomly,  X^,X2,...,Xn 
are  mutually  independent  random  variables.  Hence, 

n 

L  -  n  f(x. :K,0,C)  (17) 

i«l  1 

The  method  of  maximum  likelihood  chooses  as  estimates  those 
values  of  the  parameters  that  maximize  the  likelihood  L. 

The  procedure  is: 

1.  Take  the  partial  derivatives  of  the  natural 
logarithm  of  L  with  respect  to  each  parameter. 

2.  Set  these  equations  to  zero. 

3.  Solve  simultaneously  the  equations  for  the  values 
of  the  parameters. 
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The  Weibull  pdf  is 


K-l 


f(x:K,9,C)  =  [£Lx-l£J - ]  exp  [-(^)K], 


9 


9,K>0,x>C>0 


then 


(18) 


L  = 


Vna-Kn,  n  ,  f  Q-k  ^  ,  ~,K, 

K  9  [it  (x. -C)  ]  exp  [-9  I  (x.-C)  ] 

i=l  1  i=l  1 


(19) 


and 


ln(L)  =  nlnK-Knln9+(K-l)  Z  ln(x. -C) -9 'K  Z  (x.-C)K 

i*l  1  i*l  1 


(20) 


Taking  the  partial  derivatives  of  ln(L)  with  respect  to  K, 
9,  and  C  yields: 


3ln(L)  _  n  nlnQ+  £  ln(x.-c)-e'K  Z  (x. -C) kln(x. -C) =0  (21) 

i=l  1  i=l  1  1 


"3 


F-l- 


31n(L) 

39 


^  +  K9 


K-l 


n 


Z  (x.-C)*  =  0  (22) 

i-1  1 


31n 


=  (1-K)  Z  (xi“C)  _1  +  K9*^  Z  (x.-C)K’1  -  0  (23) 


n 


i*l 


i=l 


In  this  thesis  K  is  an  input  parameter  so  the  last 
two  equations  (22  and  23)  are  solved  simultaneously  to  deter¬ 
mine  9  and  C.  Solving  Eq  (22)  for  9  yields 


n 


K 


9  -  fj£l(xi~C)  i 1/K 
n  J 


(24) 
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The  use  of  plotting  positions  to  determine  the 
critical  values  is  based  on  the  "bootstrap"  method  (8). 


Three  plotting  positions  are  considered.  The  first, 


(i  -  .5) 
n 


(28) 


is  the  middle  of  the  interval  from  (i-l)/n  to  i/n.  The 
second  plotting  position  considered, 


i  -  .3 
n  +  .  4 


(29) 


is  the  median  ranks  plotting  position.  The  last  one  con 
sidered, 


i  i-1 
ETT  +  rTT 
- 2 - 


(30) 


is  the  average  of  the  mean  plotting  position,  i/(n+l)  and  the 
mode  plotting  position,  (i-l)/(n-l).  One  important  property 
shared  by  all  three  of  the  plotting  positions  is  the  property 
of  symmetry.  For  each  plotting  position  over  the  interval 
0  to  1 , 


i  - 


1 

1 


n 


1  - 


(n-1)  -  J- 
n 


(31) 


for  1  <  i  <  n  for  any  sample  size  n.  Johnston,  in  1980,  used 

the  median  ranks  plotting  position  in  a  method  similar  to  the 

one  used  in  this  thesis  (14:22).  At  5000  repetitions  (n  * 

5000)  there  is  no  difference  in  the  third  significant  digit 

2  2 

in  calculating  the  A  and  W  statistics  for  the  three  plotting 
positions.  For  simplicity,  plotting  position  described  by 
Eq  (28)  is  used. 

The  critical  values  are  calculated  using  linear 
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2  2 

interpolation  of  the  ordered  A  (W  )  statistics  as  the  abscissa 
axis  and  the  associated  ordinate  axis  of  ordered  plotting 
positions.  Figure  6  depicts  this  procedure.  The  arrays 
making  up  the  abscissa  and  ordinate  axes  are  composed  of 
5002  entries.  The  5000  values  of  the  plotting  position,  Eq 
(28)  with  n  =  5000  and  i  =  1 ,2  , . . .  ,5000  ,  are  entered  in  the 
ordinate  array  from  position  2  to  5001.  The  interval  [0,1) 
is  completed  by  entering  zero  and  one  into  position  1  and  5002, 
respectively.  The  5000  ordered  A  (W  )  statistics  are  entered 
in  the  abscissa  array  from  position  2  to  5001.  The  first  and 
last  entries  of  the  array  are  calculated  by  linear  extrapola¬ 
tion.  The  first  entry  is  calculated  by  linearly  extrapolating 
from  the  second  and  third  entry  (first  and  second  order  statis¬ 
tic)  subject  to  a  non -negativity  restriction.  The  5002nd 
entry  in  the  abscissa  array  is  calculated  by  linearly  extra¬ 
polating  from  the  5000th  and  5001st  entry  (second  to  last  and 
last  order  statistic) ,  not  subject  to  a  maximum  value.  The 
80th,  85th,  90th,  95th,  and  99th  percentiles  are  calculated 
by  linearly  interpolating  between  the  order  statistics  and 
the  respective  plotting  positions.  For  example,  the  80th 
percentile  is  calculated  as  shown  in  Fig.  6.  The  plotting 
position  just  greater  than  .80,  say  the  i  entry  in  the  array, 
and  the  plotting  position  just  smaller  than  .80,  the  (i-l)st 
entry  in  the  array,  are  found.  Thus,  the  80th  percentile  is 
found  by  interpolating  between  the  it^1  entry  in  the  ordinate 
and  abscissa  array  and  the  (i-l)st  entry. 
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Power  Comparison 

In  this  study  the  powers  of  the  Anderson-Darling  test 
and  the  Cramer-von  Mises  test  are  compared  with  the  modified 
Kolmogorov-Smirnov  test  and  the  Chi-Square  test.  The  power 
of  each  test  is  compared  for  several  alternative  distributions. 
The  null  hypothesis  is 


V 

Sample  deviates  follow  a  Weibull  distribution, 
shape  parameter  K  versus 

H  : 
a 

Sample  deviates  follow  some  other  distribution 

In  the  power 

study,  two  null  hypothesis  K  values  are  used. 

K  -  1.0  and 

3.5.  The  random  deviates  of  the  alternative  dis 

tributions  are  generated  using  CDC  6600  IMSL  subroutines. 

The  following  distributions  are  used  in  the  study: 

1. 

Weibull,  shape  equal  1.0 

2. 

Weibull,  shape  equal  2.0 

3. 

Weibull,  shape  equal  3.5 

4. 

Gamma,  shape  equal  1.0 

5. 

Gamma,  shape  equal  2.0 

6. 

Uniform  (1,2) 

7. 

Normal  (10,1) 

8. 

Beta  (p=l,  q«l) 

9. 

Beta  (pa2,  q»2) 

10. 

Beta  (p-2,  q*3) 

The 

Gamma  distribution  with  shape  parameter  equal 

one  is  equivalent  to  the  Weibull  distribution  with  shape  equal 
one  since  both  distributions  reduce  to  the  exponential.  The 
Gamma  with  shape  parameter  equal  two  is  presented  in  Fig.  7. 
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.5 


Fig.  7.  Gamma  Shape  *  2 

The  Beta  pdf 

f(x)  *  rEplfe]  ,  0  <  x  <  1  (32) 

for  the  above  values  of  p  and  q  are  presented  in  Figs.  8a,  8b, 

and  8c.  The  other  alternate  distributions  are  well  known  and 

will  not  be  discussed  here. 

Five  thousand  random  samples  of  size  n  are  generated 

for  each  of  the  ten  alternate  distributions.  The  respective 

2  2  2 

test  statistic,  A  ,  W  ,  K-S,  and  x  are  calculated  under  the 

null  hypothesis  that  the  random  deviates  follow  the  Weibull 

distribution  with  shape  parameter  K.  The  calculated  statistics 

are  then  compared  to  the  respective  critical  value  for  each 

2  2 

goodness-of 'fit  test.  For  the  A  and  W  tests,  the  critical 
values  generated  in  this  thesis  are  used.  The  K-S  critical 
values  presented  by  Cortes  (4:48-57)  were  generated  with  C  *  0 
and  0-1.  Thus,  the  critical  values  are  in  error  since  the 
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the  maximum  likelihood  estimators  for  C  and  9  are  not  invari¬ 
ant.  The  K-S  critical  values  used  in  this  power  study  are 
regenerated  using  invariant  estimators  of  C  and  6  and  a  simu¬ 
lation  sample  size  of  5000.  The  Chi-Square  critical  values 
used  in  this  thesis  are  generated  using  a  Monte  Carlo  simula¬ 
tion  with  five  observations  per  cell,  sample  size  n  «  25,  and 
5000  repetitions.  The  number  of  times  each  test  statistic 
exceeds  the  respective  critical  value  is  counted.  This  is 
the  number  of  times  the  null  hypothesis  is  rejected.  The 
power  of  the  test  is  the  number  of  times  the  null  hypothesis 
is  rejected  divided  by  the  total  number  of  tests,  5000.  The 
two  power  studies,  one  for  each  null  hypothesis,  are  performed 
with  sample  size  n  =  25.  The  entire  power  study  for  the  two 
null  hypotheses  and  the  10  alternate  distributions  is  repeated 
for  n  equal  to  5  and  15.  However,  for  n  equal  to  5  and  15, 
only  the  power  comparisons  between  the  A  ,  W  ,  and  K-S 
goodness-of-fit  tests  are  considered. 

Analysis  of  Critical  Values 
vs  Shape  Parameters 

The  tabled  critical  values  for  the  Anderson-Darling 
and  the  Cramer-von  Mises  tests  are  generated  for  each  value 
of  shape  parameter  K  equal  to  .5,  1.0,  1.5,  2.0,  2.5,  3.0, 

3.5,  and  4.0.  The  apparent  relationship  between  the  shape 
parameters  and  the  critical  values  is  investigated  using  re¬ 
gression  analysis.  A  graphical  relationship  between  the  shape 
parameters  and  the  critical  values  is  also  presented. 
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Computer  Programs 

The  computer  programs  used  in  this  thesis  are  pre¬ 
sented  in  Appendix  E. 
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IV.  Use  of  the  Tables 


This  chapter  discusses  the  use  of  the  tabled  criti- 
2  2 

cal  values  for  the  A  and  W  goodness-of- fit  tests  generated 
in  this  thesis  for  the  three-parameter  Weibull  distribution. 

An  example  follows  an  explanation  of  the  basic  procedures  in 
the  use  of  the  tables. 

For  both  the  Anderson-Darling  and  Cramer-von  Mises 

goodness-of-fit  tests,  a  theoretical  distribution  F(x)  is 

2  2 

compared  to  an  observed  distribution  Sn(x) .  The  W  or  A 
statistic  is  calculated  using  the  appropriate  equation  (Eq  4 
or  6).  If  the  value  of  the  statistic  exceeds  a  certain  level 
(critical  value) ,  the  theoretical  distribution  is  rejected. 

The  steps  in  applying  this  procedure  are: 

1.  Determine  the  sample  size  n  and  the  level  of 
significance,  a.  The  a-level  is  the  probability  of  rejecting 
the  null  hypothesis  that  the  sample  follows  the  theoretical 
distribution  when,  in  fact,  the  null  hypothesis  is  true. 

2.  Specify  the  shape  parameter  for  the  hypothesized 

theoretical  Weibull  distribution.  The  invariant  property 
2  2 

of  the  A  and  W  test  allows  the  scale  and  location  parameters 
to  be  unspecified. 

3.  Select,  in  a  random  manner,  the  n  observations 
from  the  population  to  be  tested.  The  n  observations  are 
ordered  from  smallest  to  largest. 
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4.  Estimate  the  unknown  location  and  scale  para¬ 
meters  from  the  observed  sample  using  the  method  of  maximum 
likelihood.  If  the  parameter  estimation  method  developed  by 
Harter  and  Moore  is  used,  care  must  be  taken  to  prevent  the 
location  estimators  from  being  mapped  onto  zero.  Any  number 
can  be  added  to  all  the  sample  values  so  the  location  esti¬ 
mator  remains  positive.  This  transformation  will  not  adversely 
affect  the  test  because  of  the  invariant  property  of  the  EDF 
statistics . 

5.  Completely  specify  the  hypothesized  Weibull 
distribution  F(x)  using  the  estimated  scale  and  location 
parameters  and  the  set  shape  parameter. 

2 

6.  Determine  the  value  of  the  test  statistic,  W 

2 

or  A  by  using  Eq  (4)  or  (6). 

7.  With  the  sample  size  n,  the  significance  level 
a  ,  and  the  shape  parameter  K,  extract  the  critical  value 
from  the  tables  in  Appendix  A  or  B. 

8.  Reject  the  null  hypothesis  if  the  value  of  the 
test  statistic  exceeds  the  critical  value.  If  the  test  sta¬ 
tistic  does  not  exceed  the  critical  value,  we  fail  to  reject 
the  null  hypothesis  and  conclude  there  is  insufficient  evidence 
to  say  the  observed  sample  does  not  follow  the  hypothesized 
Weibull  distribution  with  shape  parameter  K. 

Example 

The  following  example  illustrates  this  procedure  for 
the  Cramer- von  Mises  test. 
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The  following  ten  failure  times  are  tested  to  deter¬ 
mine  if  they  follow  a  Weibull  distribution  with  K  =  3.5  at  a 
equal  .05:  2.22,  2.26,  2.26,  2.44,  2.47,  2.49,  2.61,  2.67, 

3.14,  and  3.44  months.  The  hypothesis  tested  is: 

H  :  Sample  data  follows  a  Weibull  distribution, 

K  =  3.5  versus 

Hfl:  Sample  data  follows  some  other  distribution. 

The  Harter  and  Moore  subroutine  yields  the  estimated  para- 

A  /V 

meter  values:  location  C  =  1.427  and  scale  8  =  1.328.  Using 
these  two  estimators  and  the  shape  K  =  3.5  in  Eq  (9) ,  the 
value  of  the  hypothesized  distribution  is  calculated  for  each 
sample  value.  These  calculations  are  shown  in  Table  I. 


TABLE  I 


Example:  x  and  F(x) 


— 

i 

X 

F(x) 

i 

2.22 

.148 

2 

2.26 

.174 

3 

2.26 

.177 

4 

2.44 

.322 

5 

2.47 

.347 

6 

2.49 

.  370 

7 

2.61 

.481 

8 

2.67 

.  543 

9 

3.14 

.911 

_ “ 

3.44 

.986 

Using  Eq  (4),  the  W  statistic  equals  .144.  The  critical  value 
extracted  from  Table  XVI  with  a  equal  .05,  n  *  10,  and  K  *  3.5 
is  .142.  Since  .144  >  .142,  the  null  hypothesis  is  rejected. 


30 


Thus,  the  ten  failure  times  follow  a  distribution  other  than 
the  Weibull  distribution  with  shape  parameter  3.5. 
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V.  Discussion  of  the  Results 


This  chapter  discusses  the  results  obtained  in  this 
thesis.  The  results  as  they  pertain  to  the  four  objectives 
set  forth  in  Chapter  I  are  presented. 

Presentation  of  the  Cramer-von  Mises 
and  the  Anderson-Darling  Tables  of 
Critical  Values' 

2  2 

The  tabled  critical  values  for  the  W  and  A  goodness  - 
of-fit  tests  for  sample  sizes  n  equal  to  5,  10,  IS,  20,  25, 
and  30  and  shape  parameter  K  equal  to  .5,  1.0,  1.5,  2.0,  2.5, 
3.0,  3.5,  and  4.0  are  presented  in  Appendices  A  and  B,  res¬ 
pectively.  Each  table  presented  is  valid  for  a  specific  shape 
parameter. 

In  general,  the  critical  values  for  the  Cramer-von 

Mises  test  tend  to  monotonically  increase  for  each  level  of 

significance  as  the  sample  size  increases  from  5  to  20.  When 

the  sample  size  is  increased  from  20  to  30,  this  increasing 

trend  breaks  down.  The  trend  for  the  critical  values  for 

the  Anderson-Darling  test  is  interesting.  With  shape  para- 

2 

meter  K  equal  .5  and  1.0,  the  A  critical  values  are  monotoni¬ 
cally  decreasing  for  each  level  of  significance  as  the  sample 
size  increases  from  5  to  25.  However,  with  the  shape  para¬ 
meter  equal  to  1.5  and  greater,  the  opposite  trend  develops; 

the  critical  values  are  monotonically  increasing.  Similar  to 

2  2 
the  W  critical  values,  the  trends  for  the  A  critical  values 
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breaks  down  as  the  sample  size  becomes  large.  For  both 
goodness-of- f it  tests,  the  critical  values  may,  in  fact,  be 
approaching  a  limit  as  the  sample  size  becomes  larger  than 
25.  The  Monte  Carlo  variability  of  the  experiment  may  also 
cause  the  critical  values  to  fluctuate.  The  critical  values 
for  sample  size  greater  than  25  can  be  investigated  in  two 
ways.  First,  the  simulation  sample  size  can  be  increased 
from  5000  to  10000.  It  is  known  that  the  error  associated 
with  a  Monte  Carlo  simulation  is  proportional  to  1//TT,  where 
N  is  the  number  of  repetitions  of  the  simulation  (22:259). 

Thus,  increasing  N  would  have  the  effect  of  reducing  the  vari- 

2  2 

ability.  Second,  the  critical  values  for  both  the  A  and  W 
goodness-of-fit  tests  can  be  investigated  for  sample  sizes  n 
equal  40,  50,  and  60.  The  monotonic  decreasing  or  increasing 
trend,  thus,  may  be  easier  to  detect.  However,  either  in¬ 
creasing  the  simulation  repetitions  N  or  increasing  the 
sample  size  n  would  require  much  more  computer  time  and, 
therefore,  is  beyond  the  scope  of  this  thesis. 

Validation  of  Computer 
Programs' 

The  computer  programs  used  are  validated  by  compar- 

2  2 

ing  a  modification  of  the  A  and  W  critical  values  generated 
in  this  thesis  to  the  approximate  critical  values  calculated 
by  M. A.  Stephens  for  the  exponential  distribution  (25:359). 
This  comparison  is  possible  since  the  Weibull  distribution 
with  shape  parameter  equal  to  one  and  location  parameter 
equal  to  zero  reduces  to  the  exponential  distribution. 


Setting  K  *  1  and  C  *  0 ,  Weibull  deviates  are  gener¬ 
ated.  The  scale  parameter  e  is  estimated  and  the  critical 
values  for  the  Anderson-Darling  and  the  Cramer-von  Mises 

tests  are  calculated  for  sample  sizes  n  equal  to  5,  10,  15, 

2 

20,  and  30.  The  A  critical  values  are  modified  using  Eq  (33) 
2 

and  the  W  critical  values  are  modified  using  Eq  (34): 


A2(l  *  ^ 

-  5 1 
~~v 

(33) 

n 

W2(l  +  ■ 

.16-. 
n  J 

(34) 

7  2 

where  the  A  and  W  are  the  Anderson-Darling  and  Cramer-von 
Mises  critical  values  generated  in  this  thesis  and  n  is  the 
sample  size.  The  modified  critical  values  are  compared  to 
the  critical  values  derived  by  Stephens  for  a  level  of  signi¬ 
ficance  equal  to  .15,  ,10,  .05,  and  .01.  The  comparisons 
7  2 

for  the  W  and  A  critical  values  are  presented  in  Tables  II 
and  III,  respectively.  The  comparisons  are  good. 

Power  Investigation 

A  power  comparison  is  made  between  the  Chi-Square, 
Kolmogorov-Smirnov,  Cramer-von  Mises,  and  Anderson-Darling 
goodness-of- fit  tests  for  the  three-parameter  Weibull  distri¬ 
bution  with  the  scale  and  location  parameters  unspecified. 

The  power  comparisons  are  made  using  random  deviates  generated 
from  the  ten  alternate  distributions  discussed  in  Chapter  III. 
The  Weibull  distribution  with  shape  parameter  K  »  1  and  shape 
parameter  K  -  3.5  are  both  considered  as  the  null  hypothesis. 
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TABLE  II 


i 


2 

Cramer-von  Mises  W 


1  -  a 

W2  C 1-* 

■.  16/n) 

Stephen's 

Critical 

Values 

n 

5 

•MBM  IlMKBEt 

20 

30 

.85 

.155 

.150 

.149 

.146 

.149 

.90 

.180 

.178 

.173 

.171 

.171 

.95 

.226 

.227 

.220 

.214 

.  224 

.99 

.324 

.323 

.358 

.326 

.337 

TABLE  III 

2 

Anderson-Darling  A 


1  -  a 

A2d  +  l.S/n-5/n2) 

Stephen's 

Critical 

Values 

n 

5 

H 

.85 

.928 

.999 

.998 

.940 

.922 

.90 

1.071 

1.173 

1.153 

1.083 

1.078 

.95 

1.345 

1.492 

1.438 

1.333 

1.341 

.99 

1.996 

2.220 

2.053 

2.000 

1.957 

The  power  of  the  Chi-Square  test  is  only  included  at  sample 
size  n  *  25.  As  expected,  the  power  of  the  Chi-Square  test 

is  lower  than  the  other  three  goodness- of-fit  tests.  The 

2  2 

power  comparisons  between  the  K-S,  W  ,  and  the  A  tests  are 
discussed  by  first  considering  the  null  hypothesis  with  K  ■  1 


Shape  Parameter  K  Equal  1.0.  The  hypothesis  tested 
is 

H  :  The  sample  data  follow  a  Weibull  distribution 
0  with  shape  equal  one  versus 

Ha:  The  sample  data  follow  some  other  distribution. 

The  power  comparisons  between  the  goodness-of -fit 
tests  are  presented  in  Table  IV  for  a- level  equal  .05,  and 
Table  V  for  a-level  equal  .01.  The  power  of  the  goodness-of- 
fit  tests  when  the  null  hypothesis  is  true  should  achieve  the 
claimed  level  of  significance.  There  are  differences  for  the 
power  of  the  tests  in  the  third  decimal  place  from  the  exact 
significance  level.  These  discrepancies  can  possibly  be 
explained  by  the  Monte  Carlo  variability  of  the  experiment. 

At  the  a  equal  to  the  .05  level  of  significance,  the 

power  of  the  Cramer-von  Mises  test  is  generally  greater  than 

the  power  of  the  Anderson-Darling  or  the  Kolmogorov-Smirnov 

tests.  The  trend  is  evident  for  sample  sizes  n  equal  to  25 

and  15  and  for  the  alternate  distributions  with  the  exception 

of  the  Gamma,  shape  equal  two  when  n  ■  15.  For  the  deviates 

generated  from  the  alternate  distributions  at  ot  equal  .05  and 

sample  size  n  equal  to  25  and  15,  the  power  study  indicates 

2 

the  power  of  the  three  tests  in  decreasing  order  are:  W  , 

A  ,  and  K-S.  For  sample  size  n  *  5,  the  power  of  all  the 
goodness-of- fit  tests  are  quite  low;  they  are  all  very  close 
to  the  level  of  significance  of  the  test.  This  indicates 
the  three  goodness-of- fit  tests  are  not  very  good  when  only 
five  observations  are  used  in  the  test.  With  five 
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K-S  .012  .105  .315 


TABLE  V,  continued 


observations,  the  three  tests  are  not  able  to  differentiate 
the  alternate  distributions  from  the  Weibull  distribution 
with  shape  parameter  equal  one. 

Shape  Parameter  K  Equal  to  3.5.  The  hypothesis  tested 
is : 

H  :  The  sample  data  follow  a  Weibull  distribution 
u  with  shape  parameter  equal  3.5  versus 

Ha:  The  sample  data  follow  some  other  distribution. 

The  power  comparisons  between  the  goodness-of-fit 
tests  are  presented  in  Tables  VI  for  ct -level  equal  .05  and 
VII  for  a-level  equal  .01.  When  the  null  hypothesis  is  true, 
both  tests  do  achieve  the  claimed  level  of  significance.  The 
differences,  again,  may  be  due  to  the  Monte-Carlo  variability. 
Generally,  the  A  test  is  more  powerful  than  the  W  and  K-S 
tests.  This  study  showed  for  the  alternate  distributions 
tested  at  a  equal  .05  and  sample  size  n  equal  to  25,  15,  and 
5  the  power  of  the  three  tests  in  decreasing  order  are:  A  , 

W  ,  and  K-S.  In  contrast  with  the  first  hypothesis  (Weibull, 
shape  equal  one),  the  three  goodness-of-fit  tests  are  better 
when  the  sample  size  is  five.  However,  the  value  of  the  three 
goodness-of-fit  tests  for  small  sample  sizes  is  questionable. 

The  problem  with  small  sample  sizes  could  be  further  investi¬ 
gated  by  expanding  the  power  study  to  include  other  hypothe¬ 
sized  Weibull  distributions,  including  sample  size  n  *  10, 
investigate  additional  alternative  distributions,  or  investigate 
other  goodness-of-fit  tests.  However,  all  these  alternatives 
take  additional  computer  time  and  are  beyond  the  scope  of 
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TABLE  VII,  continued 


this  thesis. 
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Relationship  Between  the  Critical 
Values  and  the  Shape  Parameters 

The  relationship  between  the  Anderson-Darling  and 

the  Cramer-von  Mises  critical  values  and  the  Weibull  shape 

parameters  is  investigated  using  both  graphical  and  regres- 

sion  techniques.  The  graphs  of  the  W  critical  values  versus 

the  Weibull  shape  parameters  are  presented  in  Appendix  C. 

2 

The  graphs  for  the  A  critical  values  versus  the  shape  para¬ 
meters  are  presented  in  Appendix  D. 

For  the  Anderson-Darling  goodness-of-fit  test,  the 
relationship  between  critical  values  and  shape  parameters 
one  and  greater  is  investigated  using  regression  analysis. 

The  best  model  found  to  represent  the  relationship  is: 

Y  =  aQ  +  a1(l/x)  (35) 

where  Y  is  the  critical  value  and  x  is  the  shape  parameter 

greater  than  or  equal  to  one.  For  the  figures  in  Appendix 
2 

D,  the  A  critical  values  increase  from  shape  .5  to  1.0,  then 
decrease  as  the  shape  parameter  is  increased  to  4.0.  There- 
fore,  the  relationship  between  the  A  critical  values  and 
the  shape  parameter  is  only  investigated  for  shape  parameters 
1.0  and  greater.  In  Eq  (35),  ag  is  a  constant  and  a-^  is  the 
coefficient  of  the  independent  variable.  The  coefficients, 
aQ  and  a^,  are  presented  in  Table  VIII.  Separate  coefficients 
are  presented  for  sample  sizes  5,  10,  15,  20,  25,  and  30  and 
for  levels  of  significance  .20,  .15,  .10,  .05,  and  .01. 
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Coefficients  and  R  Values  for  the  Relationship  Between  the 
Anderson-Darling  Critical  Values  and  the 


2 

The  R  value  on  the  tables  is  a  measurement  of  the 
percent  of  total  variation  explained  by  the  regression  line. 

All  R  values  are  above  .80.  Twenty- five  regressions  have  an 

2  2 

R  above  .90,  and  16  have  an  R  of  .99.  Thus,  Eq  (35)  is  a 

2 

good  approximation  of  the  relationship  between  the  A 
critical  values  and  the  Weibull  shape  parameters. 

Eq  (35)  is  also  a  very  good  approximation  of  the 
relationship  between  the  Cramer-von  Mises  critical  values  and 
the  Weibull  shape  parameters  equal  to  1.5  and  greater.  For 
most  of  the  figures  in  Appendix  C,  the  critical  values  de¬ 
crease  (shape  .5  to  1.0),  then  increase  (shape  1.0  to  1.5), 
and  then  decrease  again  (shape  1.5  to  4.0).  Therefore,  the 
relationship  between  the  W  critical  values  and  the  shape 

parameters  is  investigated  for  shape  parameters  equal  to  1.5 

2 

to  4.0.  All  30  regressions  have  an  R  greater  than  .90,  and 

2 

25  are  above  .98.  The  coefficients  and  R  values  are  pre¬ 
sented  in  Table  IX, 

An  example  will  demonstrate  the  use  of  these  tables. 
Suppose  a  W  critical  value  is  desired  for  shape  parameter 
K  *  2.8,  sample  size  n  *  20,  and  an  a- level  of  .10.  The  value 
of  the  a^  and  a1  coefficients  are  extracted  from  Table  IX. 

The  coefficients  and  the  inverse  of  2.8  are  substituted  into 
Eq  ( 35) .  Thus , 

Y  *  .09  +  .099(1/2.8)  (37) 

yields  a  critical  value  of  .125. 
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Coefficients  and  R  Values  for  the  Relationship  Between  the 
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TABLE  IX,  continued 


VI.  Conclusions  and  Recommendations 


The  following  conclusions  are  made  based  on  the 
results  obtained  in  this  thesis. 

1.  The  tabled  Anderson-Darling  and  Cramer-von  Mises 
critical  values  for  the  three-parameter  Weibull  distribution 
are  valid.  In  the  Monte  Carlo  simulation,  both  tests  achieved 
the  claimed  level  of  significance  when  the  null  hypothesis  is 
true. 

2.  The  conclusions  based  on  the  power  comparison 
study  is  applicable  for  at  least  the  ten  alternate  distribu¬ 
tions  discussed  in  Chapter  III.  The  Kolmogorov- Smirnov, 
Anderson-Darling,  and  the  Cramer-von  Mises  goodness -of- fit 
tests  are  not  very  powerful  when  the  sample  size  is  five. 

When  the  hypothesized  distribution  is  the  Weibull  with  shape 

parameter  equal  1.0,  the  power  of  the  tests  in  descending 
2  2 

order  are:  W  ,  A  ,  and  K-S.  When  the  hypothesized  distribu¬ 
tion  is  the  Weibull  with  shape  parameter  equal  3.S,  the  power 

2  2 

of  the  tests  in  descending  order  are:  A  ,  W  ,  and  K-S. 

3.  The  Anderson-Darling  critical  values  are  related 
to  the  inverse  of  the  shape  parameter  for  the  Weibull  distri¬ 
bution  when  the  shape  parameter  is  one  and  greater.  The  same 
inverse  relationship  exists  for  the  Cramer-von  Mises  critical 
values  for  the  Weibull  distribution  when  the  shape  parameter 
is  1.5  and  greater. 
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Recommendations 


Based  on  observations  made  during  the  investigation, 
the  following  recommendations  are  proposed  for  further  study. 

1.  Calculate  the  critical  values  for  the  Anderson- 
Darling  and  Cramer-von  Mises  tests  for  the  three-parameter 
Gamma  distribution  when  the  scale  and  location  parameters  are 
unspecified. 

2.  Modify  Mann,  Scheuer,  and  Fertig's  S  statistic 
to  calculate  the  critical  values  for  the  three-parameter 
Weibull  distribution.  The  power  of  the  S  test  can  be  compared 
to  the  power  of  the  A  ,  W  ,  and  K-S  goodness-of-fit  tests  to 

determine  which  is  most  powerful. 

2  2 

3.  Investigate  the  A  and  W  and  possibly  other 
goodness-of-fit  tests  for  small  sample  sizes.  The  simulation 
sample  size  can  be  increased,  other  alternative  distributions 
investigated,  or  additional  shape  parameters  of  the  Weibull 

distribution  can  be  considered  as  the  hypothesized  distribution. 

2  2 

4.  Investigate  the  A  and  W  goodness-of-fit  tests 
for  large  sample  sizes.  The  critical  values  for  n  equal  to 
40,  50,  and  60  can  be  investigated  and  the  simulation  sample 
size  can  be  increased. 

5.  Investigate  the  possibility  of  developing  a  more 
efficient  (in  terms  of  computer  execution  time)  parameter 
estimation  routine  for  the  Weibull  distribution.  The  new 
routine  can  then  be  compared  to  the  Harter  and  Moore  routine 
used  in  this  thesis  to  determine  the  most  efficient. 
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Lgbles  °f  the  Cramer- von  Mises  Critical 
Values  for  the  Weibull  Distribution 
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TABLE  X 


2 

Cramer-von  Mises  W 
Shape  Parameter  Equals  .5 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.109 

.121 

.138 

.169 

.257 

10 

.118 

.133 

.157 

.195 

.293 

IS 

.120 

.138 

.160 

.204 

.310 

20 

.127 

.144 

.169 

.210 

.318 

25 

.128 

.146 

.171 

.219 

.332 

30 

.125 

.142 

.166 

.  211 

.334 

TABLE 

XI 

Cramer-von 
Shape  Parameter 

Mises  W2 
Equals  1. 

0 

■  1  1 

Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.096 

.106 

.118 

.141 

.209 

10 

.111 

.125 

.146 

.181 

.271 

15 

.114 

.130 

.150 

.186 

.279 

20 

.116 

.133 

.158 

.198 

.304 

25 

.118 

.136 

.160 

.196 

.299 

30 

.122 

.138 

.163 

.206 

.315 
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TABLE  XII 

2 

Cramer-von  Mises  W 
Shape  Parameter  Equals  1.5 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.110 

.124 

.141 

.171 

.245 

10 

.116 

.131 

.150 

.186 

.279 

15 

.117 

.133 

.153 

.193 

.289 

20 

.119 

.133 

.157 

.203 

.301 

25 

.117 

.132 

.155 

.190 

.292 

30 

.123 

.138 

.159 

.205 

.301 

TABLE  XIII 

2 

Cramer-von  Mises  W 
Shape  Parameter  Equals  2.0 


Sample 

Ci  7  A 

Level 

of  Significance 

n 

.20 

.15 

.10 

.05 

.01 

5 

.101 

.111 

.125 

.148 

.198 

10 

.104 

.118 

.134 

.165 

.  232 

15 

.104 

.117 

.136 

.163 

.237 

20 

.108 

.122 

.139 

.172 

.251 

25 

.104 

.118 

.137 

.167 

.245 

30 

.107 

.122 

.141 

.178 

.248 

TABLE  XIV 


*1 

M 


2 

Cramer-von  Mxses  W 
Shape  Parameter  Equals  2.5 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.094 

.104 

.117 

.138 

.196 

10 

.099 

.109 

.125 

.153 

.214 

15 

.098 

.110 

.127 

.153 

.203 

20 

.100 

.112 

.128 

.156 

.219 

25 

.098 

.110 

.127 

.157 

.234 

30 

.102 

.115 

.132 

.163 

.229 

TABLE 

XV 

? 

Cramer-von  Mises  W 
Shape  Parameter  Equals  3. 

0 

Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.091 

.100 

.111 

.132 

.177 

10 

.093 

.103 

.117 

.144 

.  207 

15 

.094 

.104 

.119 

.147 

.216 

20 

.095 

.106 

.124 

.150 

.212 

25 

.095 

.106 

.123 

.149 

.220 

30 

.096 

.107 

.124 

.148 

.208 

60 

i 

< 


TABLE  XVI 


Sample 

Size 

n 


Cramer-von  Mises  W 
Shape  Parameter  Equals  3.5 


Level  of  Significance 

.20 

.15 

.10 

.05 

.01 

.088 

.097 

.110 

.126 

.173 

.091 

.102 

.117 

.142 

.198 

.090 

.100 

.115 

.141 

.197 

.092 

.102 

.117 

.142 

.195 

.092 

.103 

.117 

.146 

.207 

.094 

.  107 

.124 

.153 

.218 

TABLE  XVII 


Sample 

Size 

n 


Cramer-von  Mises  W 
Shape  Parameter  Equals  4.0 


Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

.087 

,096 

.107 

.128 

.178 

.089 

.099 

.113 

.134 

.185 

.089 

.099 

.113 

.138 

.196 

.092 

.102 

.117 

.143 

.203 

.090 

.099 

.114 

.140 

.189 

.091 

.101 

.116 

.140 

.206 

APPENDIX  B 

Tables  of  the  Anderson-Darling  Critical 
Values  for  the  Weibull  Distribution 
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TABLE  XVIII 

2 

Anderson-Darling  A 
Shape  Parameter  Equals  . 5 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

1.057 

1.144 

1.309 

1.617 

2.427 

10 

.875 

.972 

1.104 

1.365 

2.090 

15 

.834 

.937 

1.075 

1.292 

1.946 

20 

.812 

.903 

1.051 

1.290 

1.974 

25 

.807 

.  888 

1.022 

1.312 

1.951 

30 

.802 

.903 

1.043 

1.304 

1.888 

TABLE  XIX 

2 

Anderson-Darling  A 
Shape  Parameter  Equals  1.0 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

1.918 

1.973 

2.040 

1.147 

2.400 

10 

1.131 

1.194 

1.284 

1.451 

1.855 

15 

.954 

1.035 

1.147 

1.355 

1.879 

20 

.874 

.961 

1.096 

1.322 

1.885 

25 

.850 

.941 

1.069 

1.287 

1.875 

30 

.838 

.927 

1.055 

1.240 

1.806 

TABLE  XX 

2 

Anderson-Darling  A 
Shape  Parameter  Equals  1.5 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.  01 

5 

.670 

.728 

.814 

.966 

1.267 

10 

.698 

.775 

.873 

1.062 

1.529 

15 

.700 

.779 

.886 

1.059 

1.482 

20 

.722 

.804 

.914 

1.131 

1.525 

25 

.741 

.822 

.939 

1.136 

1.605 

30 

.  -'31 

.807 

.923 

1.111 

1.672 

TABLE  XXI 

2 

Anderson-Darling  A 
Shape  Parameter  Equals  2.0 


Sample 

Size 

n 


Level 

of  Significance 

.20 

.15 

.10 

.  05 

.01 

.581 

.631 

.707 

.819 

1.067 

.606 

.664 

.750 

.885 

1.195 

.630 

.696 

.777 

.938 

1.279 

.635 

.698 

.791 

.957 

1.292 

.641 

.706 

.801 

.964 

1.375 

.642 

.707 

.797 

.964 

1.314 

TABLE  XXII 


2 

Anderson-Darling  A 


Shape 

Parameter 

Equals  2. 

5 

Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.542 

.588 

.647 

.  760 

.982 

10 

.573 

.623 

.698 

.819 

1.089 

15 

.581 

.645 

.719 

.867 

1.155 

20 

.594 

.656 

.742 

.885 

1.227 

25 

.601 

.657 

.743 

.899 

1.224 

30 

.590 

.655 

.734 

.881 

1.193 

TABLE  XXIII 

2 

Anderson-Darling  A 
Shape  Parameter  Equals  3. 

0 

Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.534 

.580 

.647 

.739 

.963 

10 

.558 

.612 

.683 

.804 

1.101 

15 

.  558 

.614 

.696 

.824 

1.114 

20 

.561 

.621 

.713 

.841 

1.156 

25 

.573 

.625 

.712 

.844 

1.148 

30 

.573 

.634 

.708 

.837 

1.127 
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TABLE  XXIV 


2 

Anderson-Darling  A 


Shape 

Parameter 

Equals 

3.5 

Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.515 

.561 

.621 

.714 

.904 

10 

.538 

.586 

.660 

.775 

1.044 

15 

.  542 

.595 

.672 

.809 

1.126 

20 

.554 

.602 

.683 

.  815 

1.118 

25 

.547 

.601 

.684 

.  817 

1.104 

30 

.560 

.615 

.696 

.847 

1.207 

TABLE  XXV 

2 

Anderson-Darling  A 
Shape  Parameter  Equals  4.0 


Sample 

Size 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

5 

.520 

.563 

.616 

.712 

.925 

10 

.528 

.574 

.652 

.765 

1.054 

15 

.527 

.577 

.644 

.768 

1.048 

20 

.536 

.588 

.672 

.805 

1.086 

25 

.535 

.  587 

.659 

.782 

1.099 

30 

.546 

.600 

.670 

.792 

1.122 

c.  c. 
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APPENDIX  C 

Graphs  of  the  Cramer-von  Mises  Critical 
values  Versus  the  Weibull 


Shape  Parameters 


CRITICAL  VALUES 
cP.87  0.91  0.9S  0.99 


CRIT 


2 

Fig.  10.  Shape  vs  W  Critical  Values,  Level*. 20,  n*10 
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SHAPE  VS  CRITICAL  VALUES 
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SHAPE  PARAMETER  VALUES 


Fig.  11.  Shape  vs  Critical  Values,  Level*. 20,  n=15 
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oP.88-  0.96  3.04  3.32  3.20  3.28  3.36 
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13.  Shape  vs  W2  Critical  Values,  Level*. 20,  n*25 
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SHAPE  VS  CRITICAL  VALUtS 
WEI  BULL  W2 

LEVEL=. 20  N=30 


1.00  2.00  3.00  *.00 

SHAPE  PARAMETER  VALUES 

Fig.  14.  Shape  vs  W2  Critical  Values,  Level-. 20,  n»30 


SHAPE  VS  CRITICAL  VALUES 
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Fig.  15.  Shape  vs  W2  Critical  Values,  Level”. 15,  n»5 
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Fig.  18.  Shape  vs  VT  Critical  Values,  Level*. 15, 


n*20 
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Fig.  19.  Shape  vs  W  Critical  Values,  Level*. 15,  n*25 
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Fig.  20.  Shape  vs  W2  Critical  Values,  Level-. 15,  n-30 
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Fig.  21.  Shape  vs  W  Critical  Values,  Level*. 10,  n*S 
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Fig.  22.  Shape  vs  Critical  Values,  Level=.10,  n=10 
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Fig.  24.  Shape  vs  Critical  Values,  Level=.10,  n=20 
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Fig.  25.  Shape  vs  Critical  Values,  Level*. 10,  n=25 
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.  26.  Shape  vs  W  Critical  Values,  Level*. 10 
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Fig.  27.  Shape  vs  Vr  Critical  Values,  Level3. 05,  n-5 
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Fig.  29.  Shape  vs  W2  Critical  Values,  Level*. 05,  n*15 
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Fig.  31.  Shape  vs  Wz  Critical  Values,  Level-. 05,  n«25 
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Fig.  32.  Shape  vs  W2  Critical  Values,  Level®. 05, 
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Fig.  35.  Shape  vs  W  Critical  Values,  Level®. 01,  n®15 
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36.  Shape  vs  Critical  Values,  Level*. 01,  n=20 
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Fig.  37.  Shape  vs  Critical  Values,  Level*. 01,  n=25 
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fig.  39.  Shape  vs  A*  Critical  Values,  Level-. 20,  n-5 
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Fig.  41.  Shape  vs  A2  Critical  Values,  Level*. 20,  n-15 
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Fig.  42.  Shape  vs  Critical  Values,  Level-. 20,  n=20 
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Fig.  44.  Shape  vs  A*  Critical  Values,  Level-. 20,  n-30 


Fig.  45.  Shape  vs  Critical  Values,  Level-. 15,  n-5 
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Fig.  48.  Shape  vs  A  Critical  Values,  Level*. 15,  n*20 
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Fig.  49.  Shape  vs  Critical  Values,  Level*. 15,  n=2S 
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Fig.  50.  Shape  vs  A  Critical  Values,  Level*. 15,  n=30 
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Fig.  51.  Shape  vs  A2  Critical  Values,  Level*. 10,  n-5 
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Fig.  52.  Shape  vs  A  Critical  Values,  Level*. 10,  n*10 
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CRITICAL  VALUES 


Fig.  54.  Shape  vs  AZ  Critical  Values,  Level-. 10,  n-20 
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Fig.  55.  Shape  vs  A*  Critical  Values,  Level*. 10,  n*25 
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Fig.  57.  Shape  vs  A  Critical  Values,  Level*. 05,  n*5 
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Fig.  58.  Shape  vs  A2  Critical  Values,  Level-. 05,  n=10 
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Fig.  59.  Shape  vs  A2  Critical  Values,  Level*. 05,  n-15 
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Fig.  61.  Shape  vs  A  Critical  Values,  Level*. 05,  n*25 
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Fig.  62.  Shape  vs  A  Critical  Values,  Level*. 05,  n=30 
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Fig.  64.  Shape  vs  A2  Critical  Values,  Level-. 01,  n=10 
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Fig.  65.  Shape  vs  A  Critical  Values,  Level*. 01,  n*15 


t 

i» 


125 


Fig.  66.  Shape  vs  a  Critical  Values,  Level*. 01,  n*20 
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Fig.  67.  Shape  vs  A2  Critical  Values,  Level*. 01,  n*25 
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Fig.  68.  Shape  vs  A2  Critical  Values,  Level-. 01,  n-30 
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Program  to  Calculate  the  Cramer-von  Mises  Critical  Values 


J 1 O , T2000 , I 070 , CM 1 00000 . TBOOA 1 3 , BUSH, 4022 
ATTACH, IMSL. I D=L I BRARY , SN-ASD. 

LIBRARY, IMSL. 

FTN5, ANSI-0. 

LGO. 

tEOR 

PROGRAM  CVM 

C  ********* tft*tt«**i***«t*sttttt* ***t*******.*tttMt$tttttt*t**t*t* 

C  ****************************  Mt*t**t*t.*ttttttt****t*ttt9M*tt**ttit 

C  tTHIS  PROGRAM  GENERATES  THE  CVM  STATISTICS 

C  *5000  REPS 

C  *THE  TABLES  GENERATED  ARE  VALID  FOR  THE  WE I BULL  DISTRIBUTION 
C  *N»SAMPLE  SIZE  »  5, (5) , 30 

C  *SS1«0  IF  SCALE  PARAMETER  (THETA)  IS  KNOWN 

C  *SS<«1  IF  THETA  IS  TO  BE  ESTIMATED 

C  *SS2«0  IF  SHAPE (  K  )  IS  KNOWN  ‘ 

C  *SS2«1  IF  <  K  )  IS  TO  BE  ESTIMATED 

C  *SS3»0  IF  LOCATION (C)  IS  KNOWN  j 

C  *SS3"1  IF  C  IS  TO  BE  ESTIMATED  j 

C  *C1-IN1TTAL  ESTIMATION  OF  C  (OR  KNOWN  VALUE)  1 

C  *T1-INIT:AL  ESTIMATION  OF  THETA  (OR  KNOWN  VALUE)  ] 

*EK1« INITIAL  ESTIMATION  OF  ALPHA  (OR  KNOWN  VALUE)  , 

**************************  *(.(:**% 
tttft  *%**%*%%%*%%%%*%■*%  tttttntkt***  kin.**  **xxiit*t*xt*tt*t**t't** 
COMMON / RAY / Z  < 1 00  > ,N 

COMMON /SAND /SSI , SS2, SS3, M, Cl , Tl , EKl , MR 
DOUBLE  PRECISION  DSEED 

DIMENSION  FX<60) , AA <5000) , XX (5002) , YY (5002) 

INTEGER  KEP.PP 
DSEED- I 0300. ODO 
MR-0 
REP-5002 
NOS-REP-2 
NUM-REP-2 

*  CALCULATES  <I-.5)/N  RANK 3 

YY  < 1 ) -0 
YY (REP) —1 

.  DO  405  L-2.REP-1 
YY (L) »  <  <L~I > -.5) /NOS 
405  CONTINUE 

READ*, Sal , SS2, SS3, Cl . Tl , EK1 
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PR I NT  * , SS 1 , SS2 , SS3 «  C 1 , T 1 ,  EK 1 
PRINT* 

PRINT* 

PRINT* 

PRINT  « <2X.A,F5. 1) ’ , ’SHAPE  «  ’,EKI 

PRINT  ’  <2X ,  A)  ’  ,  - - ’ 

PRINT* 

DO  100  PP*»5,30,3 

N-PP 

M-N 

DO  99  KK- 1,3000 

CALL  BGWIB<DS£ED,EK1,N,Z> 

DO  719  IK*1,N 
Zi:K>  -  1.  *2(IK)-*-2. 

719  CONTINUE 

CALL  VSRTA(Z,N) 

CALL  WE I BULL  <  CS J , TSJ , EKS  J 1 
DO  SO  L»1,N 

FX (L) =1 . -EXP <- < <Z <L> -CSJ) /TS3> **£KSJ > 

88  CONTINUE 
WCVM  »  0. 

XN-N 

DO  500  I»t,N 

WCVM  ■  WCVM  ♦  <FX < I) -<2. *XI-1 • )  t  (2.*XN>>**2 
300  CONTINUE 

WCVM  •  WCVM+1 . / < 12. *XN) 

AA<KK>  -  WCVM 
99  CONTINUE 

CALL  VSRTA  <  AA . 5000 ) 

DO  400  Lot, PEP-2  • 

XX(L+1)  ■  AA  <L> 

400  CONTINUE 

CALL  ENDPT  <XX, YY,REP,NUM> 


C 

c 

c 


*  PRINTS  PERCENTILES 

PRINT  M2X,A,  12)  ’FOR  N  -  ’,PP 

PRINT  ’  (2X.A)  ’  ,  - - ’ 

PRINT* 

DO  410  J*  30,93,3 
DO  420  II-  t.REP 
I«  REP  ♦  l-II 

IF  <YY<I> .LT. (J/100.0) ) THEN 

SLOPE  «  < YY ( 1+1 )  -  YY<IJ>  /  (XX(I+l> 

ZZ-  -SLOPE  *  XX (I >  ♦  YY< I ) 


XX<I)  ) 
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PRINT  * <2X,A, I2.A.F9.4) ’ . 

C*  THE* , J . ’ TH  PERCENTILE  IS", 

C< (J/100. >  — ZZ ) /SLOPE 
PRINT* 

GOTO  410 
END  IF 

420  CONTINUE 

41 0  CONTINUE 

DO  430  AK-l.REP 

K  -  REP+l-AK 

IF  <YY<K> .LT. .99)  THEN 

GOTO  999 

END  IF 

430  CONTINUE 

999  SLOPE  *  ( YY  <K+1 ) — YY (K) >  /  < XX <K+1 > -XX (K> > 

Zlm  —SLOPE* XX (K)  ♦  YY (K> 

PRINT  *  <2X,A.F9.4) *  , 

C’ THE  99TH  PERCENTILE  IS', 

C  < • 99—  Z  Z )  /  SLOPE 
PRINT* 

PRINT* 

PRINT* 

PRINT* 

100  CONTINUE 

PRINT  • <2X,F9.4,4X,F9.4,4X,F9.4) ’ ,  CSJ,TSJ,EKSJ 
END 
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Program  to  Calculate  the  Anderson-Darling  Critical  Values 


J20, T6000, I 070. CM1 00000 . T9006 1 3 , BUSH , 4022 
ATTACH, IMSL, ID-LIBRARY.SN-ASD. 

LIBRARY, IMSL. 

FTN5, ANSI-0. 

LGO. 

*EOR 

PROGRAM  COM 

C  *********** ******************* it* ******************************* 

C  ********************************* ******************************* 

C  *THIS  PROGRAM  GENERATES  THE  A-0  STATISTICS 

C  *5000  REPS 

C  tTHE  TABLES  GENERATED  ARE  VALID  FOR  THE  WE I BULL  DISTRIBUTION 

C  *N— SAMPLE  SIZE-  5. (S),30 

C  *SSl-0  IF  SCALE  PARAMETER  (THETA)  IS  KNOWN 

C  *SS1-1  IF  theta  IS  TO  BE  ESTIMATED 

C  *SS2-0  IF  SHAPE <  K  )  IS  KNOWN 

C  *3S2»1  IF  <  K  >  IS  TO  BE  ESTIMATED 

C  *SS3»0  IF  LOCATION (C)  IS  KNOWN 

C  *SS3«1  IF  C  IS  TO  BE  ESTIMATED 

C  <CI— INITIAL  ESTIMATION  OF  C  (OR  KNOWN  VALUE) 

C  *Ti— INI TIAL  ESTIMATION  OF  THETA  (OR  KNOWN  VALUE) 

C  *EK1«INITIAL  ESTIMATION  OF  ALPHA  (OR  KNOWN  VALUE) 

C.  **************************************************************** 

C  **************************************************************** 

COMMON/RAY/Z ( 100) , N 

COMMON / SAND / SS 1 , SS2 , SS3 , M , C 1 , T l , EK 1 , MR 
DOUBLE  PRECISION  DSEED 

DIMENSION  FX(AO) ,AA(5000) ,XX (5002) ,YY<S002> 

INTEGER  REP.PP 
DSEED- 1 0000. ODO 
MR— O 
NZERO-O 
NONE— 0 
REP-5002 
NUM-REP-2 
NOS-REP-2 
C 

C  * CALCULATES  I-.5/N  RANKS 

C 

YY ( 1 ) -0 
YY ( REP ) - 1 
DO  405  L-2, REP-1 
YY (L) - < (L-l ) -. 5) /NOS 
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403  CONTINUE 

READ*  , SS 1 , SS2 , SS3 . C 1 , T 1 ,  EK 1 
PRINT* , SSI , SS2, SS3. Cl , T 1 , EKl 
PRINT* 

PRINT* 

PRINT* 

PRINT  *  <2X,  A,F5. 1) ' , 'SHAPE  *  ’,EK1 

PRINT  •  (2X, A'  '  ,  ' - - - > 

PRINT* 

DO  100  PP-3,30,5 

N-PP 

M-N 

DO  99  KK- 1,5000 

CALL  GGWIB(DSEEF,EK1,N.Z> 

DO  719  IK-l.N 
Z  < IK) -1 . *Z  <  IK>  +2. 

719  CONTINUE 

CALL  VSRTA<Z,N) 

CALL  WE  I BULL (CSJ ,  T5J , EKSJ ) 

DO  86  L-1,N 

FX  <L)-1 . -EXP  <- < (Z  <L> -CSJ) /TSJ) **EKSJ> 

IF  <FX<L) .EQ.O. )  THEN 
FX<L)«FX<L>+.0001 
NZERO-NZERO+1 
END  IF 

IF  (FX(L).EO. 1.)  THEN 
FX«L)  »  FX(L)  -  .0001 
NONE-NONE+1 
END  IF 

ea  CONTINUE 
WAD  ■  0. 

XN-N 

DO  500  I=»1,N 
XI  -  I 

WAD-UAD+  ( 2 . * X I - 1 ) * ( LOG <  FX ( I ) )  ♦  LOG ( 1-FX (N+l-I ) > ) 
500  CONTINUE 

WAD  =  (-WAD/XN)  -  XN 
AA(KK)  «  WAD 
99  CONTINUE 

CALL  VSRTA  <  AA , 3000 ) 

PRINT* 

PRINT* 

PRINT*, 'NZERO-  ’ . N7ER0 
PRINT* , ’ NONE-  * , NONE 
DO  400  L  =  1, REP-2 
XX<L*t)  -  AA (L) 
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400  CONTINUE 

CALL  ENDPT ( X  X . YY , REP , NUM > 

*  PRINTS  PERCENTILES 

PRINT  ' <2X, A. 12) ’ , "FOR  N  -  ’ , PP 

PRINT  ’  <2X,  A)  ’  ,  ’ - ’ 

PRINT* 

DO  410  J-  00,95,3 
DO  420  II-  l.REP 
I-  REP  ♦  l-II 

IF  <  YY  ( I )  .  LT.  < ,1/ 100. 0)  )  THEN 

SLOPE  -  <YY<1+1)  -  YY(D)  /  <XX(I-t-l>  -  XX<I>) 
ZZ*=  -SLOPE  *  XX  (I)  +  YY(I) 

PRINT  ’ <2X,A, I2,A,F9.4) ■ , 

C’THE’,J,’TH  PERCENTILE  IS’, 

C< (J/100. )-ZZ> /SLOPE 
PRINT* 

GOTO  410 
END  IF 
420  CONTINUE 
410  CONTINUE 

DO  430  AK-l.REP 

K  »  REP+l-AK 

IF  <YY<K) ,LT. .99)  THEN 

GOTO  999 

END  IF 

430  CONTINUE 

999  SLOPE  -  <YY(K*1 )-YY(K> )  /  < XX <K+1 > -X X <K> > 

ZZ-  — SLOPEtXX (K)  f  YY  <K) 

PRINT  ’ <7X , A. F9. 4) ’ , 

C*  THE  99TH  PERCENTILE  IS’, 

C1.99-ZZ)  /  SLOPE 
PRINT* 

PRINT* 

PRINT* 

'  PRINT* 

100  CONTINUE 

PRINT  ’ <2X,F9.4,4X,F9.4,4X,F9.4)’,  ESJ,T5J,EKSJ 
END 


Power  Study 


J50, T2000, 1070, CM200000. TB20041 , BUSH, 4022 
ATTACH,  IMSt.,  I D-L I BR ARY ,  SN-ASD . 

LIBRARY. IMSL. 

FTN5.ANSI-0. 

LBO, PL* 10000. 

SEOR 

PROGRAM  POWER 

C  **************************************************************** 

C  tTHIS  PROGRAM  GENERATES  A  POWER  STUOV  BETWEEN  THE  FOLLOWING! 

C  *KS,  CRAMER  VON  MISES,  ANDERSON  DARLING,  AND  CHI-SQUARE  STAT. 

C  *5000  REPS 

C  »THIS  POWER  STUDY  IS  VALID  FOR  THE  WE I BULL  DISTRIBUTION 
C  *N  -  SAMPLE  SIZE  -  25 

C  *SS1*0  IF  SCALE  PARAMETER  THETA  IS  KNOWN 

C  *SS1-1  IF  THETA  IS  TO  EE  ESTIMATED 

C  *SS2«0  IF  SHAPE  <  K  )  IS  KNOWN 

C  *SS2-1  IF  (  K  )  IS  TO  EE  ESTIMATED 

C  *SS3»0  IF  LOCATION (C)  IS  KNOWN 

C  *SS?«1  IF  C  IS  TO  BE  ESTIMATED 

C  «Cl-INITIAL  ESTIMATION  OF  C  (OR  KNOWN  VALUE) 

C  *T1-INITAL  ESTIMATION  Cr  THETA  <0R  KNOWN  VALUE) 

C  *EKt* INITIAL  ESTIMATION  OF  K  (OR  KNOWN  VALUE) 

c  ************************************** ************************** 

c  **************  mmiutmiiiNtiMimmitiiitNitmitriim 

COMMON/RAY/Z  <1V0) .N 

COMMON/S AND/SSI , SS2, SS3, M, Cl , Tl ,  EK1 ,  MR 
DOUBLE  PRECISION  D5EED 
DIMENSION  FX (60) ,  FIX (60) 

C,FFX(10),ZZ(10) , CELL (5)  , 

CAAWCVM ( 5000 ) , AAWAD ( 5000  > , A AKS ( 5000  > , AACH I ( 5000  > 

INTEGER  PP 
DSEE0-20000.0D0 
MR-0 
RWCVM—O. 

RWKS-O. 

RWAD-O. 

RCHISC-O. 

NZERO-O 

NONE-O 

DO  710  IN-1,6 

FFX(IN)  -  0.  .2*  (IN-1) 

710  CONTINUE 

READ*, SSI , SS2, SS3, Cl , T1 , EK1 
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PRINT*, SSI, SS2.SS3, Cl, Tl,EKt 
PRINT* 

PRINT* 

PRINT* 

PRINT  ' <2X,A,F5. 1) * , 'SHAPE  -  ’,EKI 

PRINT  •  <2X,  A)  ’  ,  * - ’ 

PRINT* 

PP-25 

N-PP 

M-N 

DO  99  KK«=  1,5000 

CALL  GGW!8(DSEE0,2.0.N, Z) 

CO  719  IK  -  1,N 
2  (IK)  -  l.*Z(IK>  2. 

719  CONTINUE 

CALL  VSRTA  <  Z ,  N> 

CALL  WE I BULL  ( CS  J , TSJ , EKSJ ) 

DO  B8  L-l.N 

FX (L) -1 . -EXP  <- ( (Z (L) -CSJ) /TSJ) **EK3J> 

FIX <L)  -  FX(L> 

IF  (FIX (L). EG.  0,)  THEN 
FIX (L) -FIX (L) +.0001 
NZERO«NZERO+I 
END  IF 

IF  (FIX  <L> , EO. 1 .  >  THEN 
FIX (L)  «  FIX(L)  -  .0001 
NONE-NONE+l 
END  IF 

88  CONTINUE 

DO  711  IN  ■  1,3 
CELL (INI  -  0. 

711  CONTINUE 

DO  712  KI  «  2,3 

ZZ(KX>-( < ( LCG ( 1 . / ( 1 . -FF  X (KI) ) ) ) ** <1. /EKSJ) ) *TSJ)+CSJ 

712  CONTINUE 
WCVM  -  0. 

XN=N 

WAD  -  0. 

TQP=0. 0 
BOToO.O 
DO  500  I-l,N 
XI  -  I 
RL«I 

IF (RL/XN-FIX  < I )  . GT,  TOP) TOP-RL/XN-FIX (I) 
IF(FIX(I)-(RL-1)/XN  .GT.  BOT)  BOT-FIX ( I ) - (RL-1 > /XN 
WCVM  -  WCVM  +  <FX(l>-(2.*XI-l.)  /  <2.*XN>)**2 
WAD-WAD*  (2. *XI-1>  * (LOG (FIX (I) )  *  LOG < 1-Fl X (N*l-I > > ) 
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IF  (Z(I).LE.  ZZ<2)>  THEN 
CELL (I )  -  CELL ( l >  «■  1 
ELSE IF  <Z(I).L£.ZZ<3) >  THEN 
CELL  (2)  -  CELL (2)  +  1 
ELSE IF  (Z  (I  > . LE. ZZ (4) )  THEN 
CELL  < 3 )  -  CELL (3)  +  1 
ELSE IF  (Z(I).LE.ZZ(S>>  THEN 
CELL  < 4 )  -  CELL (4)  +  1 
ELSE 

CELL  <3)  -  CELL (5)  ♦  1 
END  IF 

300  CONTINUE 
DIP-TOP 

I F  < BOT  .BT.  DIF)DIF-BOT 
WKS  -  DIF 
AAKS  <KK)  -  DIF 

IF  (WKS.ST. .2071)  RWKS  »  RWKS  ♦  1  . 

WCVM  "  WCVM+1. / <12. *XN) 

AAWCVM ( KK)  -  WCVM 

IF (WCVM. BT. . 1959)  RWCVM  *  RWCVM  ♦  1 
WAD  -  (— WAD/XN)  -  XN 
AAWAD(KK)  -  WAD 

IF  (WAD. QT. 1.2873)  RWAD  -  RWAD  +  1 

CHISQ-  <  (CELL  ( 1 ) -5. >  **2) /5.  *  ( (CELL (2) -5. ) **2) /5. 

C(  (CELL  (3) -5.  )  **2)  /5.  «•  (  (CELL(4)  -5. )  **2>  /5.  + 

C( (CELL (S) -S. ) **2) /5. 

AACHKKK)  -  CHISQ 

IF (CHISQ. GT .  7.6  )  RCHISQ  -  RCHISQ  fl 
99  CONTINUE 

CALL  VSRT A ( AAKS , 3000 ) 

CALL  VSRT  A ( AAWCVM , 5000 ) 

CALL  VSRTA ( A AWAD , 5000  > 

CALL  VSRT  A  ( AACH  1,5000) 

PRINT* 

PRINT* 

PRINT* , ’ NZERO-  ’ , NZERO 
PRINT*, 'NONE-  ’.NONE 
PRINT* 

PRINT*,  ’SAMPLE  SIZE  -  ’ ,  PP 
PRINT* 

PRINT* 

PRINT*,  ’TOTAL  REJECTION  7.  FOR  K-S-  RWKS/3000 
PRINT* 

PRINT*,  ’TOTAL  REJECTION  X  FOR  WCVM-  ’,  RWCVM/5000 
PRINT* 

PRINT*, ’TOTAL  REJECTION  V.  FOR  WAD-  ’,  RWAD/5000 
PRINT*, ’TOTAL  REJECTION  Y.  FOR  CHISQUARE-  ’,  RCHIS0/3000 
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PRINT* 

PRINT*,  ' AAKS,  AAWCVM, AAWAD, AACHI * 

DO  765  IJ  -  1,5000 

PRINT  *  <2X,  I4,4(F9.4)  )  '  ,  I J, AAKS < I J > ,  AAWCVMUJ), 
C AACHI < I J) 

765  CONTINUE 

PRINT  ' <2X,F9.4,4X,F9.4,4X,F9.4>',  CSJ, TSJ, EKSJ 
END 


AAWAD  <  IJ) , 
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on  o  o  o  o  o 


Harter  and  Moore  V.'eibull  Parameter  Estimation  Routine 


COMPUTES  MLE  FOR  WEI BULL 

SUBROUTINE  WEIBULL  <CSJ . TSJ , £KSJ> 
t********************* ******* 
***************************** 

01 MENS I CN  T < 1 00 )  ,  C  < 550 ) , THET A < 550 ) , £K < 550 )  ,  X  < 36 >  ,  Y<35) 
COMMON/RAY / Z <100) ,N 

COMMON / 3 AND / SS 1 , SS2 , SS3 , M , C 1 , T 1 , EK 1 , MR 
REAL  EKSJ  * 

EN»N 

DO  2  I  »  1,N 

2  T < I ) -7  (I) 

C<1) -Cl 
THETA <11  -  T 1 
EK  ( 1 )  -  EK1 
IF  <M>  66,66,32 

32  CONTINUE 
EM=M 

31  ELNM-O. 

EP1R-MR 
MRP"  MR*1 

33  NM«N-M*1 

DO  34  1-NM,N 
EI-I 

34  ELNM=ELNM+ ALOG  >.  E I) 

IF  <MR)  66,33,74 

74  DO  75  I'M, MR 

EI-I 

73  ELNM— ELNM— ALOG (El) 

33  DO  30  J- 1.530 

IF  <J-I)  66,25,37 
37  JJ-J-1 

SK-O, 

SL-O. 

DO  6  I -MRP,  PI 

6  SK-SKXT  <I)-C<JJ)  )  **EK<JJ> 

IF  (SSI)  7.7,8 

7  THETA <J> -THETA <JJ> 

SOTO  9 

8  IF  (MR)  66,17,20 
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19  THETA (J > ■ ( (SK-MEN-EM) *  (T(M)-C(JJ)  >  **EK  <JJ>  > /EM) ** ( 1. /EK(JJ) ) 

GOTO  9 

20  X ( 1 ) “THET A  <  J  J ) 

LS-0 

DO  21  L»l,55 
UL-L-1 
LP«L-rl 
X<LP)=X(L) 

ZRK»< (T(MRP)-C(JJ) ) /X(L>  > **EK(JJ) 

Y(L)**-EK (JJ)  * (EM-EMP) /X (L)+EK<JJ) *SK/X  <L>  **  (EK(JJ>+1 . ) +EK(JJ) t (EN- 
1EM) * (T (M) -C(JJ) ) **EK(JJ> /X  (L) ** (£K(JJ) +1. >-EMR*EK(JJ> *ZRK*£XP(-ZRK 
2)/<X(L)*(l.-EXP(-ZRK) ) ) 

IF  (Y  (L) )  53,73,34 

53  LS«LS-1 

IF  (LS+L)  58,55,53 

54  LS-LS+1 

IF  <LS-L)  58,56,58 

55  X(LP)».5*X(L) 

SOTO  61 

56  X(LP)-1.5*X(L) 

GOTO  61 

58  IF  (Y(LUY(LUI)  e>0,73,59 

59  LL-LL-1 
GOTO  58 

60  X(LP)«=X  (L)+Y(L)  *  (X (L)-X<LL) ) / ( Y (LL> -Y (L) > 

61  IF  (ABS(X(LP)-X(L) l-l.E-4)  73,73,21 

21  CONTINUE 

73  THETA  (J)*=X  (UP) 

9  EK  <0 ) =£K ( JJ  > 

10  IF  <SS2>  12,12,11 

11  DO  17  I=MRP,M 

17  SL«SL+ALOG(T(I)-C(JJ) ) 

X  < 1 ) “EK < J) 

LS-0 

DO  51  L-1,55 
8LK-0, 

DO  18  I«MRP,M 

18  SLK*SLK+ (ALOG (T  < I ) -C  < JJ) ) -ALOG (THETA < J) ) > *  <T ( I > -C ( JJ) > t*X(L> 

LL»L-1 

LP-U^l 

X(LP)«X(L) 

ZRK"  <  <T (MRP) -C  < JJ ) ) /THETA  <J  >  > **X(L> 

Y (L)  ■ (EM-EMR) « (l./X <L> -ALOG (THETA ( J > > ) +5L-SLK/THETA ( J >  **X <L)+(EN- 
1EM) * (ALOG < THETA <J) ) -ALOG <T (M) -C ( JJ ) ) ) *  <T (M) -C < JJ ) ) ** X  <L> /THETA (J) 
2**X(L)+EMR*ZRK* (ALOG ( ZRK) /X  <L) >  *EXP ( -ZRK) / ( 1 . -EXP (-ZRK) ) 

IF  (Y(L>>  43,52,44 
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43  LS=LS-1 

IF  <L3+U  47,43,47 

44  LS=LS-i-l 

IF  (LS-L)  47,46,47 
43  X  (LP>  «*.  5*X  (L) 

GOTO  SO 

46  X(LP)<«1.5*X<L) 

GOTO  SO 

47  IF  (Y  <L) *Y (LL) )  49,52,48 

48  LL=LL-1 
GOTO  47 

49  X  <LP)«*X  (L)  +Y  <L)*<X(L>— X  (LL)  )  /  (Y  (LL)  —  Y  (L)  ) 

50  IF  <AI?S<X<LP)-X<L>)-l.E-4>  52,52,51 

51  CONTINUE 

32  EK<J)«XU.P> 

12  C<J)*C<JJ) 

62  IF  <S53)  23,25,14 

14  IF  < 1 . -EK  <  J )  )  1  A, 78, 78 

78  IF  (SS1+SS2)  57,57,16 

16  X(1)*=C<J> 

LS=0 

DO  23  L-1.55 
SKl-O. 

SR«=0. 

DO  13  I-MRP.M 

SK1«SK1+<T<I)-X (L) ) *» (EK  ( J) -1. ) 

13  SR-SR+1./<T<I)-X'L>> 

LL-L-1 

LP-L*1 
X<LP)»X (L) 

2RK**  < ( T  < MRP ) - X <L> > /THETA <J>  >**EK(J> 

Y (L) “ < 1 • — EK <  J ) ) 4SR+EK<J) * (SK1+ <EN-EM> *  <T CM> -X  <L>  >  ** (EK(J)-l. ) > 

1 /THETA <J)**EM.J)-EMfi*EK<OU*ZRK*EXP<-ZRK>/  l  <7  (MRP)-X  <L)  )  *  (l.-EXP 
2<-2RK> ) > 

IF  <Y<U)  39,24,40 

39  LS-LS-l 

IF  (LS+L)  70,41,70 

40  LS=LS+1 

IF  <LS-L>  70,42,70 

41  X<UP)».5*X <L) 

GOTO  22 

42  X  <LP)  *>.  3*X  <L)+.5*T <  1 ) 

GOTO  22 

70  IF  <Y<L>*YO.L>)  72,24,71 

71  LL-LL-t 
GOTO  70 
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I 


72  x(lp>»x<l)-<-v(L)*<x<u-x;ll>>/(V<li.)-y(l)) 

22  IF<ABS'X<LP>-X(L))-:.E-4>  2',, 24, 23 

23  CONTINUE 

24  C(J>*»X(LP> 

GOTO  25 

57  C(J>«T(1) 

25  IF (MR)  66.38,69 

38  DO  63  >i,M 

1F(C<J)t1. E-4-T ( I ) )  68,67,67 

67  MR-MR+1 

63  C ( 1 >  «T ( 1) 

68  IF (MR)  66,69,31 

69  SK-O. 

SL-O. 

DO  36  I=>MRP,M 
SK»SK+(T(I>-C(J>  >**EK(J> 

36  SL-SL+ALO!3<T<n-C(J>  ) 

ZRK-( (T(MRP)-C(J) > /THETA <J)) t*EK(J) 

EL-ELNM+ ( EM-EMR ) * ( ALOG (EK ( J  > ) -EK ( J ) * ALOG'c  THETA ( J  > )) ♦ <EK  M  > - 1 . > *SL- 
1 (SK-*-  (EN-EM) *  <T (M)-C(J) >  *»EK<0) ) / (THETA (J ) *TEK (J >  > +EMR*ALOQ ( 1 . -EXP 
2<-ZRKl ) 

IF ( J-3)  30,27,27 

27  IF (ABS (C(J)-C(JJ) ) -1 . E-4)  28,28,30 

28  IF(ABS(THETA(0)-THETA(JJ) )-l.E-4)  29,29.30 

29  IF (ASS (EK (J)-EK (JJ> ) -1 . E-4)  4,4,30 

30  CONTINUE 

4  CONTINUE 

CSJ»C(J) 

TBJ-THETA(J) 

EKSJ-EK(J) 

66  RETURN 

END 

♦EOR 

I*  ,0a  ,0,  ,0.  ,  X,  ,  1. 

♦EOR 

♦EOF 


Program  to  Evaluate  the  Endpoints 


* SUBROUTINE  TO  EVALUATE  ENDPOINTS 

SUBROUTINE  ENDPT (XX, YY, REP, NUM> 

INTEGER  REP 

DIMENSION  XX  <REP) ,  YV!REP> 

SL0PE«<YY(2)-YY<3) I  /  <XX(2)-XX(3>> 

B  «  YY (2)  -  SL0P£*XX<2) 

VI  *  -B/SLQPE 

IF  (Vl.LT.O.)  THEN 

Vl-O. 

END  IF 
XX(1)«V1 

SLOFE*  < YY  < NUM )  -YY  < NUM+ 1 )  )  /  <XX  (NL'M) -XX  (NUM+’. )  > 
Bl-YY(NUM)  -  SLOPE* XX !NUM) 

V2-U.-S1  ) /SLOPE 

XX<REF)-V2 

RETURN 

END 
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